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ABSTRACT 

We present computer simulations of interferometric imaging with the VLT interferometer and the AMBER instru- 
ment. These simulations include both the astrophysical modelling of a stellar object by radiative transfer calculations 
and the simulation of light propagation from the object to the detector (through atmosphere, telescopes, and the 
AMBER instrument), simulation of photon noise and detector read-out noise, and finally data processing of the 
inter ferograms. The results show the dependence of the visibility error bars on the following observational pa- 
rameters: different seeing during the observation of object and reference star (Fried parameters ro, bject=2.4m, 
^o,ref.=2.5m), different residual tip-tilt error (£ tt ob j ect =2% of the Airy disk diameter, £ tt ref =0.1%), and object 
brightness (if object =3. 5 mag and 11 mag, if re f.=3.5 mag). Exemplarily, we focus on stars in late stages of stellar 
evolution and study one of its key objects, the dusty supergiant IRC +10420 that is rapidly evolving on human 
timescales. We show computer simulations of VLTI interferometry of IRC +10 420 with two ATs (wide-field mode, 
i.e. without fiber optics spatial filters) and discuss whether the visibility accuracy is sufficient to distinguish between 
different theoretical model predictions. 
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1. INTRODUCTION 
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The Very Large Telescope Interferometeiu (VLTI) with its four 8.2 m unit telescopes (UTs) and three 1.8 m auxiliary 
telecopes (ATs) will certainly establish a new era of optical and infrared interferometric imaging within the next 
Qh" few years. With a maximum baseline of up to more than 200 m, the VLTI will allow the study of astrophysical key 
objects with unprecendented resolution opening up new vistas to a better understanding of their physics. 

The near-infrared focal plane instrument of the VLTI, the Astronomical MultiBEam RecombinerBi (AMBER), 
will operate between 1 and 2.5 /im and for up to three beams allowing the measurement of closure phases. In a second 
phase its wavelength coverage is planned to be extended to 0.6 /im. Objects as faint as K = 20 mag are expected to 
be observable with AMBER when a bright reference star is available, and as faint as K = 14 mag otherwise. 

Among the astrophysical key issuesS are, for instance, young stellar objects, active galactic nuclei and stars in 
late phases of stellar evolution. The simulation of a stellar object consists in principle of two components: 

(i) the calculation of an astrophysical model of the object, typically based on radiative transfer calculations 
predicting, e.g., its intensity distribution. To obtain a robust and non-ambiguous model, it is of particular 
importance to take diverse observational constraints into account, for instance the spectral energy distribution 
and visibilities. 

(ii) the determination of the interferometer's response to this intensity signal, i.e. the simulation of light propagation 
in the atmosphere and the interferometer. 

Often, only one of the above parts is considered in full detail. The aim of this study is to combine both efforts and to 
present a computer simulation of the VLTI performance for observations of one object class, the dusty supergiants. 
For this purpose, we calculated a detailed radiative transfer model for one of its most outstanding representatives, 
the supergiant IRC +10 420, and carried out computer simulations of VLTI visibility measurements. The goal is to 
estimate how accurate visibilities can be measured with the VLTI in this particular but not untypical case, to discuss 
if the accuracy is sufficient to distinguish between different theoretical model predictions and to study on what the 
accuracy is dependent. 



2. THE SUPERGIANT IRC +10 420: EVOLUTION ON HUMAN TIMESCALES 



The star IRC +10 420 (= V 1302 Aql = IRASJ.9244+1115) is a^putstanding object for the study of stellar evolution. 
Its spectral type changed from F8I+ in 1973B to mid-A todayHl3 corresponding to an effective temperature increase 
of 1000-2000 K within only 25 yr. It is one of the brightest IRAS objects 4u£-tp its very strong infrared excess by 
circumstellar dust and one of the warmest stellar OH maser sources knoimBBilJ Large mass-loss rates, typically of 
the order of several lO _4 M0/yr were determined by CO observations .tilH IRC +10 420 is believed to be a massive 
lumincms^staTj .currently being observed in its rapid transition from the red supergiant stage to the Wolf-Rayet 
phase .13113'tHlBLl Its massive nature (initially ~ 20 to 40 M ) can be concluded from its distance (d = 3-5 kpc) and 
large wind velocity (40 km/s) ruling out alternative scenarios. IRC +10 420 is the only object observed until now in 
its transition lathe Wolf-Rayet phase. The structure of the circumstellar environment of IRC +10420 appears to be 
very complex J13_and scenarios proposed to explain the observed spectral features of IRC +10 420 include a_rotating 
equatorial disk,li3 bipolar outflows,E3 and the infall of circumstellar material onto the star's photosphere. 113 

Several infrared speckle and coranographic observationsB00@0 were conducted to study the dust shell of 
IRC +10420. The most recent studyEil reports the first diffraction-limited 73 mas bispectrum speckle interferometry 
of IRC +10420 and presents the first radiative transfer calculations that model both the spectral energy distribution 
and the visibility of this key object. In the following we will briefly describe the main results and conclusions of this 
study which will serve as astrophysical input for the VLTI computer simulations presented in the next section. 

2.1. Visibility measurements and bispectrum speckle interferometry 

Speckle interferograms of IRC +10420 were obtained with the Russian 6 m telescope at the Special Astrophysical 
Observatory on June 13 and 14, 1998.EH The speckle data were recorded with our NICMOS-3 speckle camera 
through an interference filter with a centre wavelength of 2.11 jiva and a bandwidth of 0.19/im. The observational 
parameters were as follows: exposure time/frame 50 ms; number of frames 8400; 2.11 jiva seeing (FWHM) ^l'/O. A 
diffraction-limited image of IRC +10420 with 73 ma^psalution was reconstructed from the speckle interferograms 
using the bispectrum speckle interferometry method.E3'E§Ed The modulus of the object Fourier transform (visibility) 
was determined with the speckle interferometry methodES Figure [l] (left panel) shows the reconstructed 2.11 /im 
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Figure 1. Left: Azimuthally averaged 2.11 /im visibility of IRC +10 420 with error bars for selected frequencies. 
This visibility function consists of a constant plateau (visibility ~ 0.6) caused by the unresolved central object and 
a triangle-shaped low- frequency function caused by the faint extended nebula. Right: Azimuthally averaged radial 
plots of the reconstructed diffraction-limited 2. 11 /mi- images of IRC +10420 (solid line) and HIP 95447 (dashed line). 

visibility function of IRC +10420. There is only marginal evidence for an elliptical visibility shape (position angle 
of the long axis ~ 130° ± 20°, axis ratio ~ 1.0 to 1.1). The visibility 0.6 at frequencies > 4cycles/arcsec shows that 
the stellar contribution to the total flux is ~60% and the dust shell contribution is ~40%. The GauB fit FWHM 
diameter of the dust shell was determined to be ^fwhm = (219 ± 30) mas. By comparison, previous 3.8 m telescope 



K-band observations^ found a dust-shell flux contribution of ~50% and ^fwhm = 216 mas. However, as will be 
shown later, a ring-like intensity distribution appears to be much better suited than the assumption of a Gaussian 



distribution whose corresponding FWHM diameter fit may give misleading sizes (see Sect. 2.2). The right panel 
of Fig. |l] displays the azimuthally averaged diffraction-limited images of IRC +10420 and the unresolved star HIP 
95447. The if -band visibility will strongly constrain radiative transfer calculations as shown in the next section. 

2.2. Dust-shell models 

The spectral energy distrubion (SED) of IRC +10420 with 9.7 and 18 /im silicate emission features is shown in Fig. ||. 



It corresponds to the '1992' data set used by Oudmaijer et alfl) and combines VRL (October 1991), near-infrared 
(March and ApriL1992) and Kuiper Airborne Observatory (June 199JJ photometry^ with the IRAS measurements 
and 1.3 mm data.ES Additionally, we included data for A < 0.55 /im.E3 In contrast to the near-infrared, the optical 
magnitudes have remained constant during the last twenty years within a tolerance of ~ 0™1. 

IRC +10 420 is highly reddened due to an extinction of Ay tal w 7 m by the interstellar medium aiicLthe circum- 
stellar shell. From polarization studies an interstellar extinction of Ay & 6 m to 7 m was estimated.EIlEj Based on 
the strength of the diffuse interstellar bapds, E(B — V) = 1™4 ± 0™5 can be inferred for the interstellar contribution 
compared to a total of E(B — V) =j=j 2™4E2 We will use an interstellar Ay of 5 m adopting the method of Savage & 



MathisES with A v = 3AE(B - V)\ 




Figure 2. SED (top) and 2.11 /im visibility (bottom) for a superwind model with Y = rjr\ = 4.5 and different 
amplitudes. The inner shell obeys a 1/r 2 density distribution, the outer shell p-l/r 1,7 density distribution. Model 
parameters are: black body, T e ff = 7000 K, T\ = 1000 K, ro.55 Mm = 7.0, silicates,l23 standard grain size distribution^!! 
with a max = 0.45 /im, and F Q ut = 10 4 . The symbols refer to the observations (see text) corrected for interstellar 
extinction of A v = 5 m . 




Figure 3. Left: Fractional contributions of the emerging stellar radiation as well as of the scattered radiation and 
of the dust emission to the total flux as a function of the wavelength for a superwind model with Y = rjr\ = 4.5, 
A = 40 and a 1/r 1,7 density distribution in the outer shell. Model parameters are: black body, T e ff = 7000 K, 
Ti = 1000 K, 70.55/Ltm = 7.0, silicates ,ti2l standard grain size distribution!^ with a max = 0.45 /im, and Y out = 10 4 . 
Right: Corresponding normalized intensity vs. angular displacement The (unresolved) central peak belongs to 
the central star. The inner hot rim of the circumstellar shell has a radius of 35 mas, and the cool component is 
located at 155 mas. Both loci correspond to local intensity maxima. 

In order to model both the observed SED and 2.11 /im visibility, we performed radiative transfer calculations for 
dust shells assuming spherical symmetry. We used the code DUSTYE3 which solves the spherical radiative transfer 
problem utilizing the self-similarity and scaling behaviour of IR emission from radiatively heated dust.L2l To solve 
the radiative transfer problem including absorption, emission and scattering several properties of the central source 
and its surrounding envelope are required, viz. (i) the spectral shape of the central source's radiation; (ii) the dust 
properties, i.e. the envelope's chemical composition and grain size distribution as well as the dust temperature at 
the inner boundary; (iii) the relative thickness of the envelope, i.e. the ratio of outer to inner shell radius, and the 
density distribution; and (iv) the total optical depth at a given reference wavelength. The code has been expanded 
for the calculation of synthetic visibilities. 

We calculated various models of the dust shell of IRC +10420 considering black bodies and model atmospheres as 
central sources of radiation, different silicates, grain-size and density distributions. We refer to Blocker et alEil for a 
full description of the model grid. The remaining fit parameters are the dust temperature, Ti, which determines the 
radius of the shell's inner boundary, 7*1, and the optical depth, r, at a given reference wavelength, A re f. We refer to 
A re f = 0-55 jam. Models were calculated for dust temperatures between 400 and 1000 K and optical depths between 
1 and 12. Significantly larger values for r lead to silicate features in absorption. 

The near-infrared visibility strongly constrains dust shell models since it is, e.g., a sensitive indicator of the grain 
size. Accordingly, high-resolution interferometry results provide essential ingredients for models of circumstellar 
dust-srplls. In the instance of IRC +10 420 (assuming the central star to be a black body oilT e ff = 7000 K), the 
silicatet3 grain sizes, a, were found to be in accordance with a standard distribution function,El n(a) ~a -3,5 , with 
a ranging between <2 m i n = 0.005 //in and <2 m ax — 0.45 /im. 

However, the observed dust shell properties cannot be fitted by single-shell models but seem to require multiple 
components. At a certain distance we considered an enhancement over the assumed \/r x density distribution. The 
best model for both SED and visibility was found for a dust shell with a dust temperature of 1000 K at its inner radius 
of 7*1 = 69 R*. At a distance of r = 308 R* (Y = r/r\ = 4.5) the density was enhanced by a factor of A = 40 and 
and its density exponent was changed from x = 2 to x = 1.7. These fits for SED and 2.11 jam visibility are shown in 
Fig. ^| The various flux contributions at 2.11 /im are 62.2% stellar light, 26.1% scattered radiation and 10.7% dust 
emission (see Fig. ||), i.e. the radiation emitted by the circumstellar shell itself consists of 71% scattered radiation and 
29% direct dust emission. The shell's model intensity distribution is shown in Fig. [| and was found to be ring-like. 



Table 1. Parameters of simulations 



wavelength 


brightness 


opt. throughput 


read-out noise 


quantum efficiency 


exposure 


photon number 


2.11 /im 


K = 3.5 mag 


0.108 


15 e~ 


0.6 


50 ms 


4.42 • 10 6 



This appears to be typical for optically thin shells (here ro.55 Mm = 7, T2.n Atm = 0.55) showing limb-brightened dust- 
condensation zones. Accordingly, the interpretation of the observational data by FWHM GauB diameters may give 
misleading results. The ring diameter is equal to the inner diameter of the hot shell 69 mas), and the diameter of 
the central star amounts to ~ 1 mas. The bolometric flux, Fboi, is 8.17- 10 -10 Wm -2 corresponding to a central-star 
luminosity of L/L Q = 25 462 • (d/kpc) 2 . 

This two-component model can be interpreted in terms of a termination of an enhanced mass-loss phase roughly 90 
yr (for d = 5kpc) ago. The assumption that IRC +10420 had passed through a superwind phase in its recent history 
is in line with its evolutionary status of an object in transition from the Red-Supergiant to the Wolf-Rayet phase. 
The mass- loss rates of the components can be determined to be M\ = 7.0 10 -5 M /yr and M2 = 1.1 10 -3 M /yr. 

3. INTERFEROMETRY WITH THE VLTI AND THE AMBER INSTRUMENT 

In the previous section, the spatial intensity profile of the dusty supergiant IRC +10 420 (see Fig. j|) was derived by 
means of radiative transfer models and their comparison with photometric and interferometric observations. This 
2.11 /im intensity profile will serve as object intensity profile in the simulation of monochromatic VLTI observations. 
The next steps are the simulation of light propagation from the object to the detector (through atmosphere, telescopes, 
and the AMBER wide-field mode instrument), simulation of photon noise and detector read-out noise, and finally 
data processing of the inter ferograms. A schematic view of Michelson interferometry with the VLTI and the AMBER 
camera is given in Fig. ||. 

3.1. Computer simulation of interferometric imaging 

Fig. H shows a flow chart of our simulation of interferometric imaging with the VLTI (ATs, or UTs with adaptive 
optics) and the AMBER camera in the wide-field mode (i.e. without fiber optics spatial filtering). In a first step, 
an array of Gaussian distributed random numbers is generated and convolved with the correlation function_pf the 
atmospheric refraction index variations in order to generate wavefronts degraded by atmospheric turbulence .LJ The 
typical size of the atmospheric turbulence cells is given by the Fried parameter 7*0 . After the simulation of the 
entrance pupil, the next step incorporates the tip-tilt correction but we allow for a residual tip-tilt error 5 tt . In 
the next step a typical Michelson output pupil is simulated (pupil reconfiguration). The output pupil is chosen 
such that (i) in the optical transfer function the off-axis peaks are separated from the central peak, and (ii) the 
inter ferograms are sampled with the smallest number of pixels to assure the lowest influence of detector noise. The 
following step includes light propagation through the beam combiner lens to the focal plane. The squared modulus of 
the Fourier transform of the complex amplitude in front of the beam combiner lens yields the intensity distribution of 
a Michelson interferogram of a point source. In the next step the required object intensity distribution is simulated 
(given here by the intensity distribution of IRC +10 420) to obtain the Michelson interferogram of the object: The 
Fourier transformation of the object intensity distribution, calculated at those spatial frequencies covered by the 
the simulated interferometer baseline vector, is multiplied with the off- axis peaks of the transfer function of the 
generated Michelson interferogram.Eil Finally, Poisson photon- noise and detector read-out noise is injected to the 
inter ferograms. The noise level depends, among other parameters, on the number of detectable photons, the total 
optical throughput of the interferometer and the quantum efficiency of the detector. Details^J are shown in Table |l|. 

3.2. Visibility data for IRC +10 420 

We performed simulations of IRC +10 420 AT visibility observations in the AMBER wide-field mode and studied 
the influence of various observational parameters on the visibility accuracy. Visibility error bars were, for example, 
obtained for the following observational parameters: different seeing during the observation of object and reference 
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Figure 4. Schematic view of interferometric imaging with the VLTI AMBER camera. 
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Figure 5. Flow chart of the simulation of VLTI AMBER interferograms (wide-field mode, i.e. without fiber optics 
spatial filtering). 



star (Fried parameters ro, bject=2.4m, ro, re f. = 2.5m), different residual tip-tilt error (£ tt)0b j ect =2% of the Airy disk 
diameter, ^tt,ref.=0.1%), and object brightness (if object =3. 5 mag and 11 mag, if re f.=3.5 mag). In the computer 
experiments (a)-(c), object and reference star were assumed to have the same brightness (if =3. 5 mag, see Table [l]), 
in experiment (d) fainter objects (if =11 mag) were simulated as well. Fig. || shows the results of the simulations 
(a) to (c) based on the A = 40 intensity profile (see Fig. ||) for baselines of 50 and 100 m together with the model 
predictions for IRC +10420. To obtain error bars each simulation was repeated three times. The insetted table lists 
the parameters of the simulations. 

Simulation (a) based on 2400 interferograms represents the ideal case of excellent seeing conditions (ro, bject=2.5 m, 
r*o,ref.— 2.5 m) and an almost perfect tip-tilt correction (residual tip-tilt error S tt = 0.1%). The visibility error AV 
amounts to ±0.0029. Simulation (b) illustrates the influence of a larger residual object tip-tilt error (£ tt =2%) leading 
to AV ± 0.0036. Simulation (c) shows the impact of different seeing conditions for object and reference star. In the 
last simulation (d) of this series (shown only in the insetted table of Fig. ^) an IRC -f 10420-like intensity distribution 
is assumed but the simulated if -magnitude is 11 mag (instead of 3.5 mag). Although the photon number decreases 
to N = 4.41 • 10 3 (i.e. by a factor of 1000, see Table fj]), the visibility accuracy is still high, i.e. AV ± 0.0036. 

4. CONCLUSIONS 

We have presented computer simulations of interferometric imaging with the VLT interferometer and the AMBER 
instrument in the wide- field mode. These simulations include both the astrophysical modelling of a stellar object 
by radiative transfer calculations and the simulation of light propagation from the object to the detector and sim- 
ulation of photon noise and detector read-out noise. We focussed on stars in late stages of stellar evolution and 
examplarily studied one of its most outstanding representatives, the dusty supergiant IRC +10420. The model in- 
tensity distribution of this key object, obtained by radiative transfer calculations, served as astrophyiscal input for 
the VLTI/ AMBER simulations. 

The results of these simulations show the dependence of the visibility error bar on various observational parame- 
ters. With these simulations at hand one can immediately see under which conditions the visibility data quality would 
allow us to discriminate between different model assumptions (e.g. the size of the superwind amplitude). Inspection 
of Fig. [6] shows that in all studied cases the observations will give clear preference to one particular model. Therefore, 
observations with VLTI will certainly be well suited to gain deeper insight into the physics of dusty supergiant s. 
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Figure 6. Dependence of the error bars of simulated AT- VLTI/ AMBER (wide-field mode) observations of 
IRC ±10420 at 2.11 fim on (i) seeing differences of object and reference star observations (Fried parameter, ro, differ- 
ences); (ii) different residual tip-tilt errors, 5 t t> of object and reference star observations; and (iii) object brightness. 
Lines refer to radiative transfer models of different superwind amplitudes A. The symbols refer to AT- VLTI/ AMBER 
simulations (each with 7V=2400 interferograms) of the A=40 intensity profile (see Fig. ||) for baselines of 50 and 100 m 
(100m: <7=227 arcsec -1 ). To better distinguish between the simulations, the data points belonging to one baseline 
are somewhat shifted with respect to the spatial frequency. The error bars are based on a triple iteration of each 
simulation. The table gives the parameters of the simulations (a) to (d). Note that though simulation (d) corresponds 
to an IRC ±10420-like intensity distribution, it refers to a much fainter object of if=llmag (instead of 3.5 mag). 
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